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We report on ^^As-NMR measurements in single crystalline BaFe 2 (Aso. 96 Po.o 4)2 for magnetic fields parallel to the 
orthorhombic [110]o and [100]o directions above the structural transition temperature Ts - 121 K. A large difference 
in the linewidth between the two field directions reveals in-plane anisotropy of the electric field gradient, even in the 
tetragonal phase. This provides microscopic evidence of population imbalance between As-4px and 4py orbitals, which 
reaches In^ - nylfn^ + ny\ ~ 15% at T ^ Ts and is a natural consequence of the orbital ordering of Fc-Sd^z and dy^ 
electrons. Surprisingly, this orbital polarization is found to be already static near room temperature, suggesting that it 
arises from the pinning of anisotropic orbital fluctuations by disorder. The effect is found to be stronger below ~ 160 K, 
which coincides with the appearance of nematicity in previous torque and photoemission measurements. These results 
impose strong constraints on microscopic models of the nematic state. 


The relationships among spin, orbital, and lattice degrees 
of freedom are a central issue in iron pnictides since fluc¬ 
tuations in these channels are potentially at the origin of 
high-Tc superconductivity. The problem is best illustrated 
by the presence of two interrelated structural and magnetic 
transitions, at temperatures Ts and Tn, respectively, which 
both disappear near the optimal T^ in most pnictides.Al¬ 
though the tetragonal-to-orthorhombic transition slightly pre¬ 
cedes (or coincides with) the antiferromagnetic (AFM) transi¬ 
tion, the latter is not necessarily a consequence of the former: 
whether the orthorhombic distortion is driven by magnetic or 
orbital fluctuations has been the subject of continuing debate^ 
(see Ref. 5 and references therein). Furthermore, the tetrag¬ 
onal phase above Ts turns out to be unconventional: various 
electronic properties have been found to break the fourfold 
symmetry of the Fe-As planes The relative role and the 
microscopic origin of the spin and orbital instabilities lead¬ 
ing to such an “electronic nematic state” are also controver¬ 
sial.^^’Actually, because magnetic and orbital degrees of 
freedom are entangled and most likely cooperate,^^ it is chal¬ 
lenging to determine which one, if any, is dominant. While 
both the magnetic and orbital scenarios find support in exper¬ 
iments 14,20,39^3 nematic state remains puzzling. 

Here, we report on NMR measurements above Ts in under¬ 
doped BaFe 2 (Aso. 96 Po. 04 ) 2 - We demonstrate that orbital po¬ 
larization of the As-4p orbitals, related to Fe-3J polariza¬ 
tion, is present within electronic domains in the tetragonal 
phase, even without any applied uniaxial stress. We estimate 
the magnitude of this polarization and show that it is static. 
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We further reveal an unanticipated temperature dependence: 
while the onset of the nematic state, as inferred from torque 
and photoemission experiments,^^’is manifested by an up¬ 
turn near 160 K in our linewidth data, static orbital polariza¬ 
tion is found to be already present at much higher temper¬ 
atures. This suggests that static short-range orbital order is 
first nucleated around defects and that its evolution towards 
longer-range order is directly involved in the manifestation of 
nematicity at the macroscopic scale. 

Single-crystalline BaFe 2 (Aso. 96 Po. 04)2 was synthesized by 
the conventional self-flux method.^'^ Edges of the crystal were 
cut along [110]o in the orthorhombic notation (As-Fe direc¬ 
tion), [100]o (Fe-Fe direction), and [ll0]o [see Fig. 1(a)] 
and the resulting dimensions of the crystal were 2.0 x 1.4 x 
0.25 mm^. For the following NMR measurements, the field 
was parallel to the aZ?-plane to within +1°. ^^As-NMR spec¬ 
tra were obtained by sweeping the frequency in a fixed field 
of 15 and 6 T. The structural transition atTs - 121 K 
was determined from a kink in the resistivity measured on 
a crystal from the same batch while the magnetic ordering 
at^N - 121 K was determined from the sharp drop in the 
NMR signal intensity, signifying the simultaneous structural 
and magnetic transition Ts = T^. 

Since ^^As nuclei have spin I = 3/2, the spectrum 
shows three lines corresponding to the transitions f = 
m m - 1 (m = ±112,312) dX the frequencies 

(6, 0) = (r/27r)yUo//[l + K(6, 0)] + (m - l/2)y(^, 0) + 
(2“^-order quadrupolar correction), where (6, 0) indicates po¬ 
lar and azimuthal angles of the applied field in the orthorhom¬ 
bic basis (ao,bo,c), yl27i = 7.2919 MHz/T is the ^^As gy- 
romagnetic ratio, and K(6,(p) represents the Knight shift, 
which is proportional to the magnetic susceptibility = 
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Fig. 1. (Color online) (a) VtPn layers and applied field directions. Left 
panel: the central As site without P dopant among its NN is named Asq. Mid¬ 
dle and right panels: the central As sites having one P dopant among their 
NN in the [lOOJo and [OlOJo directions are named and re¬ 

spectively. As-NMR central line (b) and satellites (c) at 200 K for different 
field orientations. Continuous lines are fits to Asq (blue) and Asi (red) sites. 
Expected resonance positions are shown as dashed lines. 



Fig. 2. (Color online) (a) T dependence of the quadrupole frequency |y|. 
(b) T dependence of the distribution Ay/|y| where Ay is the full-width at half¬ 
maximum of a satellite line. Fitting error bars are within the symbols. ^^As 
satellites at 250 (c) and 130 K (d) fox H || [lOOJo and [llOJo, respectively. 
Broken blue and red lines represent Asg and ASq sites, respectively, obtained 
from Gaussian fits to the lineshape (their sum is shown as a solid green line). 
The two sites are identical for // || [llOJo, while they split for // || [lOOJo- 
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(d) T= 130 K 
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0,oj = 0). The angular dependence of the quadrupole fre¬ 
quency Va is expressed as y(6, 0 ) = (Vcl2)(m - 1 /2)(3 cos^ 6 - 
1 - T] sin^ 6 coslc/)). is proportional to the diagonal com¬ 
ponents of the electric field gradient (EFG) tensor Vaa = 
d^VIdrl'. Va = eVaaQI2h (a = a,b,c), where Q is the elec¬ 
tric quadrupole moment. The EFG tensor of the so-called 
Aso sites far from any P dopant [see Fig. 1(a)] is diagonal¬ 
ized in the orthrhombic-basis (Gq, bo, c).^^ The in-plane EFG 
anisotropy is defined as 77 = \Vaa - Vbb\/\Vcc\ = - Vb^/Wcl 

with Va + v/? + Vc = 0. Therefore, 7 / = 0 in a local tetragonal 
environment. 

Small additional peaks observed in the tail of the central 
line [Fig. 1(b)] are attributed to those As sites (hereafter called 
Asi) having one P dopant among their four nearest neighbors 
(NNs). The Asi resonance splits into two peaks for ^ || [100]o 
because Asj^^^^° and Asj^^^^° sites [defined in Fig. 1(a)] are 
inequivalent, while they are equivalent for if || [110]o and 
[ll0]o. In Ba(Fei_jcCojc) 2 As 2 , NNs to Co dopants produce 
similar peaks.^’^^’'^^ The angle dependence of the resonance 
frequency of the Asi center peak is reproduced with the pa¬ 
rameters |y^| = 5 JMHz, j] = 0.24, and S = 30°, where S is 
the tilt angle of the principal axis (z-axis) of the EFG from the 
crystalline c-axis [Fig. 1(b)] to [100]o. The Knight shift of Asi 
is assumed to be the same as that of Asq. These parameters 
enable us to estimate the positions of the Asi satellite peaks 
[Fig. 1(c)]. As seen in the satellite peaks for || [100]o, weak 


and broad tails are observed. We consider that the broad tails 
arise from the Asq site next to Asi in the doped P and Asi 
directions, since these tails have the same angle dependence 
as the Asi site. It is clear from Fig. 1(c) that the central part 
of each satellite signal essentially arises from the Asq sites 
where the effect from the P-dopant is small. 

Below, we discuss the electronic state of the Asq sites 
by analyzing two quantities: (i) the absolute value of their 
quadrupole frequency |y|, directly determined from half of 
the separation between the high- and low-frequency satellites 
|v(0)| = 1 / 3^1 - /_i^_ 3|/2 = ly^Kl + 77 COs 20)/2, (ii) the 
distribution of quadrupole frequencies Ay/|y|, where Ay is the 
full-width at half-maximum obtained by the Gaussian fit of 
the Aso satellite peak. 

The central observation of our study is shown in Figs. 2(a) 
and 2(b). While the values of |y| are almost identical for ^ || 
[110]o and for H || [100]o, the values of Ay/|y| for these two 
orientations are different and, most importantly, this differ¬ 
ence increases on cooling. The difference in the line width is 
independent of the analysis procedure. The contrasting behav¬ 
ior between H \\ [110]o and H || [100]o is also observed in the 
line width of the ^^As-NMR center-line arising from the Asq 
site measured in 6 T, as shown in Fig. 3. We investigated the 
origin of the linewidth difference from the field dependence 
of the linewidth, and confirmed that the central linewidth is 
not determined by the hyperfine interaction, but mainly by the 
EFG. Whereas the linewidth decreases for ^ || [110]o, that the 
width for ^ || [100]o is anomalously large and that it grows 
with decreasing T can only be explained if there are two unre¬ 
solved sites for if II [100]o having identical EFG tensor com- 
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Fig. 3. (Color online) T dependence of the ^^As-NMR arising from the 
central transition obtained for// || [llOJo (a) and [100]o(b). (c) T dependence 
of the full-width at half-maximum of the centerline from the Asq site in // || 
[llOJo and [llOJo- A/center for // || [lOOJo starts to increase below 160 K, as 
observed in the linewidth of the satellite peak shown in Fig. 2 (b). 


ponents but a nonzero rj and in-plane principal axes perpen¬ 
dicular to each other, although no clear splitting was observed 
in the satellite spectra at 130 K. The experimental results can 
be understood by considering the presence of perpendicular 
domains, in each of which the in-plane fourfold symmetry of 
the EFG is broken. Note that the central and satellite lines 
were found to have identical widths for if || [110]o and for 
[lT0]o (Fig. 1 (b) and 1 (c)). 

Thus, to estimate the EFG anisotropy parameter rj, we per¬ 
formed new fits to the sharp part of the satellite peaks {i.e., 
excluding the tails since the P-dopant effect is observed) for 
H II [100]o with two Gaussian functions, each having a 
linewidth fixed to the value measured for [110]o [Figs. 2(c), 
2(d)]. The fiat-top shape of the spectra is well reproduced 
by the two- Gaussian peaks. The corresponding two sites are 
named As^ and As^ and they are defined by their respective 
quadrupole frequencies Iv^l < |y/,| for || [IOOJq. It is impor¬ 
tant to realize that the frequency difference between the satel- 



Fig. 4. (Color online) (a) T dependence of in-plane EFG anisotropy pa¬ 
rameter 77 for X = 0.04 (this work) and x = 0 (from Ref. 45). For x = 0.04, 
Tj is derived from the angular dependence of the satellite linewidth with re¬ 
spect to the field. For x = 0, however, it is deduced from the satellite splitting 
owing to twinned orthorhombic domains.If domains are present above Ts, 
this method does not allow the detection of 77 ^ 0 locally, (b) T dependence 
of Va,b for = 0-04 together with the x = 0 data from Ref. 45. is the 
quadrupole frequency estimated from band-structure calculations, is the 
value calculated using a point charge model and the lattice parameters from 
Refs. 50 and 51. (1 - 7 oo)yJfJ, for x = 0 is derived from the y vs. plot 
(see Supplemental Material). 


life peaks of Asq and Asq actually corresponds to Iv^ - vt^I for 
perpendicular domains in which the a- and b- axes are rotated 
by 90°. By dividing this value by IvJ = Iv^ + vt,! = 2|y[iio]J, rj 
can be estimated. As shown in Fig. 4(a), rj is already nonzero 
at 250 K but it increases on cooling. 

In general, the EFG has two origins: y^ = v®” + (1 - 7 oo)v|fS 
where y°” is the contribution from the on-site charge density 
and yj^^ is the contribution from the surrounding charged ions 
in the lattice, multiplied by the antishielding factor (1 -yoo)-^^ 
Therefore, in-plane anisotropy (rj ^ 0) can occur because of 
nematic order, (i.e ., electronic states locally break the fourfold 
symmetry) and/or the locally orthorhombic lattice symmetry 
(while still globally tetragonal). 

In order to obtain information on the on-site wave function, 
we examine the relative weights of these two contributions. 
We performed band-structure calculations for v = 0 within 
the local density approximation (FDA) for non-spin-polarized 
BaFe 2 As 2 ^^ and using the lattice parameters reported previ¬ 
ously.^^ Since the obtained y^^^^ = +1.95 MHz at 300 K is 
similar to the experimental value ly^,^! = 1.5 MHz,"^^ it is 
found that the on-site contribution is dominant, and the sign 
of y^^T, is chosen to be positive in Fig. 4(b) for both v = 0 
and 0.04 in the tetragonal phase. On the basis of theoretical^^ 
and experimental^^ grounds, we estimate the T -independent 
lattice contribution (1 - yoo)y^^\ = 0.45 MHz from the extrap¬ 
olation of the linear relationship between y and the spin part 
of the Knight shift (with T as an implicit parameter) 
to ^ 0 for the parent compound v = 0 (see Supple¬ 
mental Material).This contribution thus accounts for less 
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Fig. 5. (Color online) T dependence of the orbital polarization \nx - Uy \ for 
X = 0 and 0.04 estimated by and Eq. (2). The right inset shows x = 0.04 
data on an enlarged scale. The black arrow shows a crossover in the T de¬ 
pendence, coinciding with observations of a nematic state in bulk measure¬ 
ments.^^ The left inset is a schematic description of Fe-3J and As-4/7 orbital 
orderings. 


than half of the total Va,b in both x = 0 and 0.04. In addition, 
the observed nonzero r] value above Tg is much larger than rj 
estimated from the point-charge model using the lattice dis¬ 
tortion reported in Ref. 22. Taking these results into account, 
the main origin of the in-plane anisotropy is due not to local 
orthorhombic distortion, but to the twofold symmetry of the 
electronic states. 

is then deduced by subtracting (1 - yoo)y^a ffoni the 
experimental values of Va for v = 0 and 0.04. By defining 
i = (x, y, z) with 4px and 4py orbitals of As directed toward As 
NN, y®” can be expressed in terms of the occupation number 
Hi for each As-4p/ orbital,^"^’and the rii and the orbital po¬ 
larization corresponding to an occupation difference between 
As-4px and 4py are estimated from 




^(y,z,x) + ^(z,x,y)\ 


( 1 ) 


and 


\n.-ny\=—K"-v‘^, 


( 2 ) 


where Vo = = 56.90 MHz with (r ^> 4 ^ = 

6.958 a.u. Note that this simple modeling neglects any spatial 
modulation of the amplitude of the orbital polarization and 
the possibility that it is zero for some of the nuclei (in which 
case we should fit with three lines in an all-or-nothing model). 
Therefore, the value of \nx - %| should only be considered 
as an order-of-magnitude estimation. In particular, \nx - ny\ 
could be overestimated at high temperatures and thus it would 
extrapolate to zero at a significantly lower temperature than 
what the inset of Fig. 5 suggests. 

\nx-ny \ for v = 0.04 is much smaller than in the orthorhom¬ 


bic phase of V = 0, but it clearly increases with decreasing 
r, especially below ~ 160 K. This is clearly observed in 
the line width of the central transition shown in Fig. 3. This 
characteristic temperature of ~ 160 K is consistent with the 
temperature at which twofold symmetry appears in the pre¬ 
vious magnetic torque, lattice parameters, and Fe-3 J orbital 
polarization.^^’However, more work is needed to determine 
whether there is a continuous change in the magnitude and 
length scale of the orbital polarization (our data can be fit to a 
Curie Weiss law) or whether there is a sharp change at ~ 160 
K. 

While quantitative evaluation of the polarization of the Fe- 
3d orbitals has been reported in a few experiments for the 
AFM ordered state, no such measurement has been re¬ 

ported for the As-4p orbitals above Ts, to the best of our 
knowledge. Furthermore, the low time scale of NMR (set by 
the inverse linewidth of ~ 0.5 MHz) implies that the observed 
orbital ordering is static. The imbalance between 4px and 4py 
occupations most likely results from the ordering within the 
Fe-3 J orbitals, that is, the degeneracy of Fe-3 and 3dyz or¬ 
bitals is lifted. Indeed, LDA calculations revealed that a re¬ 
duction of the ¥e-3dyz partial density of states (PDOS) is ac¬ 
companied by a decrease in the As-4py PDOS in the AFM 
ordered state for v = 0.^^ This suggests that orbital ordering 
above Ts for BaFe 2 (Aso. 96 Po. 04)2 produces rix < riy for the As- 
4p orbitals. This in-plane orbital anisotropy must participate 
in lifting the degeneracy of stripe spin correlations between 
[100]o and [010]o directions. A similar relationship between 
the lattice symmetry breaking and the enhancement of spin 
fiuctuations has been observed in LaFeAsO.^^’^^ 

An important and unanticipated aspect of our data is that 
p at the Aso sites, albeit small, is nonzero even at the high¬ 
est measured temperature of 250 K, that is, well above the 
claimed nematic transition at T* ~ 160 We suspect that 
orbital ordering would be entirely fiuctuating at high temper¬ 
atures in the absence of disorder but it is made static owing to 
pinning by defects (most likely P dopants) and/or microstrain 
induced by the dopants.Such a phenomenon is analogous 
to the pinning of charge-density-wave (CDW) fiuctuations 
above the CDW phase transition^^’. In iron-based super¬ 
conductors, there is actually theoretical^^’and experimen¬ 
tal evidence from scanning tunneling spectroscopy^^’and 
transport measurements^^ that defects play a role in the ne¬ 
matic response. We stress that intrinsic spin and/or orbital ne¬ 
matic correlations are not produced by defects but they can be 
made visible by them. We expect this to be particularly true in 
the absence of any aligning field, i.e., in the unstressed tetrag¬ 
onal phase. In a simplified picture, static, but short-range cor¬ 
related, orbital order is nucleated around P sites. 

To summarize, our NMR work revealed unequal pop¬ 
ulations of the As 4px and 4py orbitals owing to Fe- 
3d orbital ordering in the (unstressed) tetragonal phase of 
BaFe 2 (Aso. 96 Po. 04 ) 2 - The data reveal that these orbital-nematic 
correlations are involved in the appearance of a nematic state 
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detected below ~160 K by other probes.Furthermore, the 
magnitude of the orbital polarization, its static nature, and its 
persistence at temperatures well above 160 K place strong 
constraints on microscopic models of nematic order in this 
class of superconductors. 
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Supplemental Material for 
^‘Emergent of Orbital Nematicity in the 
Tetragonal Phase of BaFe2(Asi_;,P;,)2’’ 

Majumder et al. reported the linear relationship between 
and y with an implicit parameter T > Ts in iron pnictides.^^ 
It was indicated from theoretical point of view that the tem¬ 
perature dependence of v can be influenced by the spin sus¬ 
ceptibility in the presence of the mode-mode coupling be¬ 
tween charge and spin density fluctuations in itinerant mag¬ 
nets.Such linear relationships are actually observed in the 
various itinerant magnet systems like MnSi,^^ PrCoAsO,^^ 
NaCo 204 ,^^ LiFeAs,^^ and BaFe 2 As 2 .^^ Since the existence 
of the mode-mode coupling between these two channels 
seems to be reasonable, the contribution of (1 - yoo)y^a 
V = 0 can be extracted by extrapolating the linear portion of 
Va versus plot to ^ 0. For the estimation, we used 
Vc and Xc^^” adopted from Ref.^^ to derive (1 - 7 oo)v|f^ (Fig. 
6). Then, (1 - which is the value used in Fig. 4(b), 

is derived by multiplying -1/2 to (1 - 7 oo)v|f^ by following 
the relationship = -(1/2) x from Laplace’s equation 
y}f^ + = 0 and yjf^ = yjf^ above Ts. Finally, we obtain 

(1 - 7oo)y^^i, to be 0.45 MHz. 
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